Quantum Dynamics in Non-equilibrium Strongly Correlated Environments 
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We consider a quantum point contact between two Luttinger liquids coupled to a mechanical 
system (oscillator). For non- vanishing bias, we find an effective oscillator temperature that depends 
on the Luttinger parameter. A generalized fluctuation-dissipation relation connects the decoherence 
and dissipation of the oscillator to the current-voltage characteristics of the device. Via a spectral 
representation, this result is generalized to arbitrary leads in a weak tunneling regime. 



The quest to build a scalable quantum computer has 
recently lead to a series of spectacular experiments where 
macroscopic quantum states were coherently manipu- 
lated and measured . These experiments for the first 
time give an opportunity to study the effects of indirect 
continuous measurement Q on an individual quantum 
system. Understanding the measurement process is not 
only important for the development of a quantum com- 
puter, but also is a fundamental problem of quantum 
mechanics. In some cases, theoretical analysis of mea- 
surement process based on explicit models describing the 
coupling between the system and electrical measurement 
apparatus has been performed 1^-^. 

Correlations play a crucial role in protecting quantum 
coherence in macroscopic systems and enabling manip- 
ulation of the quantum states. In sub- micrometer elec- 
tronic systems, the Coulomb interaction becomes impor- 
tant and can lead to Coulomb blockade, a subject of in- 
tensive research both in the contexts of "classical" and 
quantum-coherent electronic devices. The more subtle ef- 
fects of itinerant electron-electron interactions, however, 
have not been studied in the context of quantum mea- 
surement. These effects are often important since in or- 
der to interface with a quantum device, at least a part 
of the apparatus has to be scaled down to the device 
size. The study of such correlation effects is therefore 
important in developing the understanding of a realistic 
quantum measurement. 

To analyze the role of correlations in the measurement 
apparatus, we study here a specific example of two Lut- 
tinger liquid leads electrically coupled to a quantum sys- 
tem, with the tunneling current being influenced by a co- 
ordinate of the system. The Luttinger liquids are a gen- 
eralization of the Fermi liquids to one dimension, where 
electron-electron interactions lead to dramatic renormal- 
ization effects near the Fermi surface [|lO|. The interac- 
tions in the Luttinger liquids can be parametrized by the 
dimensionless constant g, which describes repulsive inter- 
actions if g < 1, attractive if <? > 1, and a non- interacting 
Fermi liquid for g = 1. One experimental realization of 
Luttinger liquids is carbon nanotubes We develop 

the formalism for an arbitrary quantum system, which 
can be, for instance, a two-level system (qubit), or a 
quantum oscillator. For g > 1, we find that the effect 



of the measurement is equivalent to coupling to a heat 
bath, with an effective temperature reduced relative to 
the Fermi liquid case The effective temperature de- 
pends on the density of states of the Luttinger liquids 
at the tmmel contact, Tcs — {a + l)^^qV/2, where the 
current / oc V"'^^. The exponent a depends on the tun- 
neling geometry and on the Luttinger parameter g. This 
provides an interesting example in which the Luttinger 
parameter determines, not an exponent, but a prefactor 
for a universal expression which is of first order in the 
bias voltage. It had been hoped that the two-terminal 
conductance of a one-dimensional wire would be propor- 
tional to g I p^ , providing another such example, but this 
turns out not to be the case . Via a spectral represen- 
tation, we find a relation valid for any leads in the weak 
tunneling regime with tunneling particles of charge q 



Tefi'-(q//2)(dF/d/). 
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Examples of experimental realizations of our model as 
applied to the measurement of a quantum oscillator are 
shown in Fig. |. In Fig. 1(a) and Fig. 1(b), the tun- 
nel junction is formed by a nanotube (Luttinger liquid) 
and a metal (Fermi liquid), while in Fig. 1(c) and 1(d) 
both sides of the junction are Luttinger liquids. In the 
first example. Fig. 1(a), the tunnel current between the 
gate and the nanotube is used to monitor the transverse 
nanotube oscillations. The characteristic oscillation fre- 
quency is about 1 GHz for a 100 nm nanotube [Q, which 
makes it possible to achieve the quantum regime at about 
50 mK. In Fig. 1(b), a short stiff nanotube in the STM 
mode [ p^ is used to perform vibrational spectroscopy of 
an adsorbate loosely bound to a metal surface. The posi- 
tion of the atom/molecule modulates the tunnel current. 
In the last two examples a kink in the nanotube formed 
either mechanically or due to a 5-7 defect plays the role 
of the tunnel contact 16|. The presence of the defect 
will lead to formation of a localized optical phonon mode 
above the nano-tube phonon band (above 200 meV p^t), 
that will couple to the tunnel current by modifying the 
tunneling matrix element. Alternatively, the chemical 
kink defect can be functionalized by adsorbing an atom 
or molecule |Q, whose vibrations will also modify the 
tunneling between the two Luttinger legs. 

Tunneling Between Two Luttinger Liquids — We con- 
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sider the problem of tunneling between two Luttinger 
liquids, when the tunneling is coupled to an external sys- 
tem, such as a quantum oscillator, or a spin. The Hamil- 
tonian is 

n = nT + nf + + na, (2) 

where Wq is the Hamiltonian for the measured system, 
referred to as an oscillator, Tif 2 the Luttinger liquid 
Hamiltonians for the two leads, with Luttinger parameter 
g and with a potential difference ^ = qV. We define the 
electron tunneling operator T{t) = = 0,t)'^2{x — 

0, t) + h.c, where ^'1^2(2^, t) are fermion operators in the 
leads. The term includes c-number terms as well as 
operators that do not commute with Tip- 

Following Kane and Fisher p6| , we consider tunnel- 
ing via a weak link, that is tunneling between two semi- 
infinite leads with x = located at the ends of the leads. 
We consider the case of repulsive interactions, g < 1- For 
spinless fermions, the end density of states of a lead at 
energy e is proportional to e"™"*, with aendC?) — l/ff— 1- 
For carbon nanotubcs, where the fermions have spin 
but the interactions arc only in the charge sector, one 
finds acndCs) = (l/ff — l)/4 Then, the exponent 

a for a tunnel junction between two leads with gi,(?2 
is a = aend(ffi) + aend(ff2)- For tunneling between 
two infinite leads, the end density of states is replaced 
with atun = (.9 + .9^^ ^ 2)/2 for spinless fermions and 
tttun = {g + g^^ — 2)/8 for nanotubes. 




(a) (b) 




(c) (d) 

FIG. 1. Experimental realizations of the model. 

We use a Keldysh formalism Our procedure 

closely follows that used to study noise in Luttinger liquid 



tunneling |2C|] . Let us suppose that initially the oscillator 
and leads are decoupled, with density matrices po for the 
oscillator and ^^2 ^or the leads, so that the full density 
matrix p — pq ® pi ® P2 ■ For now, we assume that the 
leads are at zero temperature. After the interaction be- 
tween systems is turned on at time t — —oa^ the systems 
become coupled. Define the scattering operator S by 

-i J Hdt -i J Hdt 

S = TcC °° e ~~ =1, (3) 

where the operator denotes time ordering along the 
Keldysh contours. Points along the forward branch 
(—00 — > 00) are ordered with increasing times, while 
those in the return (00 — > —00) are ordered with decreas- 
ing times, with those in the return branch ordered after 
those in the forward branch. We will occasionally use a 
superscript /, r on t to indicate to which contour t be- 
longs. The expectation value of any product of operators 
0{ti),0{t2)^ ... can be obtained by 

{0{t,)0{t2)...)=TT[p(T,0{t{)0{ti)...s)], (4) 

where we work in the Schroedinger representation 
throughout. 

Tracing out the Luttinger liquids, we obtain the new 
scattering operator 

-i J Tiodt -i J Hadt 

vSoff = Tc e - e -~ X (5) 

±i J dudt2 n(ti)o(t2) (;i°;4''-/j)''+'^ +■■■ 

e , 

where the plus sign is chosen before the integral over 
ti, t2 if ti, t2 are in different branches and the minus sign 
is chosen if they are in the same branch. The ± sign 
in the denominator of the exponential is taken positive 
if t2 is after ti and negative if t2 is before ti. We have 
used (T(ti)r(t2))L = 2cos(M(ti -i2))/(e±i(i2-ii))"+', 
where the expectation value {)l for the leads is the expec- 
tation value for decoupled Luttinger liquids at zero tem- 
perature [20I] . In this expectation value, there is an ad- 
ditional factor dependent on the density of states, which 
may be absorbed into the normalization of il. The ... 
in Eq. (|^) denote terms higher order in fi. Assuming 
ri^/i" << 1, these higher order terms may be neglected 
by a power counting, valid for a > 0; for a = 0, there 
may be logarithmic infrared divergences which upset this 
naive power counting, and in this case we must also as- 
sume a sufficiently large voltage to neglect these terms. 
The expectation value of an operator 0{t) becomes 

(0(<)) = Tr{po(rcOW5eff)} = (OW), (6) 

where we have defined O = 



2 



{0{t{))L T / dun{t){0{t{)T)i 



(7) 



Here the ... denote connected expectation values which 
are higher order in fl, and where =F is chosen negative for 
t on the forward contour and positive for t on the return 
contour. The operator O depends only on the oscillator 
coordinates and not on the leads. 

The exponential in Eq. involves a product 

n(ti)0(t2), where ti,t2 may be on either the for- 
ward or return contour. We introduce fl'^iti) — 

n{t{) + n{t{),n''{ti) = n{t{) - n{t{). Then, 



-217"(ii)r2"(t2)6l(ti -t2)Im 



2cos(^(ii - ^2)) 
(e + *|t2-ii|)"+2 
2cos(/ifti - ^2)) 



(8) 



Again assuming that f2^/i" << 1, we can make a fur- 
ther simplification in the exponential of Eq. (|^) , by using 
the Bloch-Redfield approximation |2^, that fl°''^{t2) ~ 
^mo{t2-ti)^a,s(^^-^^-z-Ho{t2-ti)_ Corrections to the Bloch- 
Redfield approximation arise if operators f2 are inserted 
between ti,t2; such corrections to S^s will be of order 
1^2 — ti\U^ n"'^^ . Let us write ft = where i,j denote 
eigenstates oi 7io with energies Eij. Then, integrating 
over t2, 



-t J Hodt -i J Hodt 

5cff = T,e ~ e — X (9) 



hi 



n"(t)n?(t)A(Ai,_Ei-Bj)-Sl°(i)Sl°(t)S(/j,£;i-_Ej) 



where we define Sifi, AE) = 2 cos {fit)e'^^*Re[{e + 
i\t\)-"'-^]dt, A{^i,AE) = -2/°^2cos(Ait)e'^-^*Im[(e-f 
dt. The terms in n'^Q" in Eq. (|) produce deco- 
herence, and are equivalent to averaging over a randomly 
fluctuating field coupled to fl, while the terms ilTi^ pro- 
duce dissipation. 

Taking \Ei — Ej\ < /i, so that the correct poles in the 
integrals for A, S are determined by the sign of fi, one 
finds that 

S{^l,AE)^^(^Ii^l + AE)+Ii^^-AE)'^ (10) 

A{^l, AE) = i (/(^ + AE) - - Ai?)) + 

where ... denotes imaginary terms, possibly singular as 
e — > 0, which may be absorbed into a renormalization 
of Ho, and hence dropped. We have defined I{fi) — 
27r^"+Vr(a + 2). Eqs. (|Jl^) are the main results. 

Average Current and Noise — Here ql{fj.) is equal to 
the current flowing at = 1. We now recom- 
pute the current within the present formalism, in or- 
der to obtain corrections to the current due to fluc- 
tuations in f2. The current operator at time t-^ = 



is JiQJ') = qifliO^ ) {^\{0, 0^)'^' 2(0,0^) - h.c.). From 
Eqs. (|p); the leading contribution to (J) is of order fi^, 
±iq JZ^ dt{n{t)n{0f))2 sin (lit) /{e±it)°'+^. This van- 
ishes when integrated over t on the forward contour, so 
we can assume that t is on the reverse contour. Ap- 
plying the same Bloch-Redfield approximation, we get 
q J^^ dte'^E,-'^^^* {n{t)n{0f))2 sin {^it) /(e + ^^)"+2. Do- 
ing this integral yields 



{j)^q{a,,io)n{o)) 



27r 



r(a-f 2) 



{fi + E,~E,)'^+\ (11) 



We now consider fluctuations in the current, 
{J{t{)J{t2)) ■ The order fl^ contribution to the current- 
current correlation function, from Eq. (0) , is given by 



,2 



q'mDQiti)) 



ie + i\t2-ti\)^+"' 



(12) 



This represents the shot noise in the tunnel junction 
slightly modulated by oscillator. To next order in fi^, 
from Eq. (^, we must compute 



-±±y J dt3dt4n{t{)n{t^^)n{t3)n{t4)) X (13) 

— 00 

{T{t{)T{ti)T{t,)T{t,))L. 

Even for c-number 57, the calculation of Eq. (|l3|) is in- 
volved in this case, the calculation yields a result of 
order /i^"r2''|t2 — plus terms with lower powers of 

fi. However, for operator f2, there is a contribution from 
Eq. ( p^ ) which is of order /z^""*"^, which will dominate 
over the previous contribution for \t2 — ti\ >> /i^^, the 
time regime we now consider. 

The expectation value {T{t{)T{t{)T{ts.)T{t4)) l = 
{T{t{)T{t:i))L{T{t{)T{tA))L +t3^t4 + connected. The 
last term, a connected expectation value of 4 T-operators, 
gives the contribution of order /i^"r2'*|i2 — to 
Eq. (p^). We ignore this, and consider only the other 
terms. Integrating over t^^t^, and applying a similar 
procedure to that used to calculate the current above, 
we arrive at the following contribution to {J{ti)J{t2)): 

lHT<^)Hf^ + E.~ i?,)"+^(M + E, - £;,)"+! (14) 

{n,,{t\)njk{tl)^{t{)n{t{)). 

For the expectation value {T(t{)T{t'i)) L{T{t{)T{t4)) l to 



be significant, ti — t^ 



t2 - ti ^ ^ , justifying 



the application of the Bloch-Rcdfield approximation. For 
\t2 — ti\ >> this approximation docs not change 

the time ordering of t^jt^. As claimed, Eq. ( [l^ ) is of 
order reflecting a modulation of the current by 

the oscillator. Eq. ( |l^ ) decays on a time scale of order 
the inverse damping coefficient, I/7, of the oscillator. In 
the case of c-number 17, this term is neglected: when 
computing (J J) — (J) (J), it cancels. 
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Density Matrix — For AE << /z, it is possible to write 
the results above in a compact form. Define p{t) to be 
the density matrix for the oscillator at given time t. Then 

use Eqs. (||J|) to compute {0{t^)) for any operator O; the 
result is be a linear differential equation for the expecta- 
tion values. Then use {0{t^)) = Tr{p{t)0) to derive the 
equation for the density matrix: 



A^S + 2I) 



1 dl r 

2d^ 



n 



,{a,p} 



/ r 

2 



(15) 



with A = [Hq,H]. Comparing to the results for Fermi 
liquid leads the effective temperature of the oscillator, 
determined by the ratio of the third (decoherence) term 
to the second (dissipation) term in Eq. (^, is 



off 



{a + l)-\V/2. 



(16) 



For the current, we find (J) = q{Vi^)I{p) + g(Ari}d//d/i. 

Spectral Representation — While these results were de- 
rived for Luttinger liquid leads, they are more general. 
Assuming sufficiently small J7, we find Eq. (^ with 
2cos(^(ii — t2))/{e ± i{t2 — ti))"^^ replaced by the ap- 
propriate expectation value in the leads, l(r{t\)T(ti))L- 
Let the density of states (particle and hole, respec- 
tively) at energy E be p^{\{E) in leads 1,2 respec- 
tively. Then, defining pP{E) = J dEi p\{Ei)pll{E - Ei), 

p'^{E) = JdEi p\{Ei)pI{E - El), the expectation value 
is /d£;(pP(£;)e=F»(^+^')(ti-t2) 4. p''(^)eT«(s-A')(ti-t2))^ 

with the minus sign chosen if ti is after t2 and the plus 
sign otherwise. Then going through the same steps, we 
find Eq. (§|,^o|) in general, with /(/i) replaced by the 
appropriate current-voltage characteristic of the device: 
I{fj,) = 2'K{pP{-n) + p'^{p)). The relation between S,A 
and / generalizes the fluctuation-dissipation relation de- 
rived between noise and current [ p4[ . For A << /i, we 
arrive at Eq. ([|). 

Applications — We now consider the specific case of 
a harmonic oscillator with frequency loq and mass m, 
linearly coupled to the tunneling, f2 = f2o + ex. We 
consider a regime for which the fig term dominates and 
the oscillator coordinate x only weakly modulates the 
current. For p, > ujo, the average current is 



j(m) 



q[ml 



c^S^ 



c^A 



2mLUoA 2mujQ 



(17) 



where / = I{p)] S — S{p, ujq); A = A{p,, cjo)- 

The noise spectrum can be evaluated 
from Eqs. ( |l^ , [l^ ). Fourier transforming Eq. ( [l^ ) gives 
the shot noise contribution 



dte*'^*(r!(o)r!(t)}L(T(o)r(i))j 



2q^^ll, (18) 



for 1^1 <C \p\- The modulation, Eq. (13), to leading order 
in c^ from Eq. (0), yields 



CJq)2 -I- "f'^Uj'^ 



(19) 



where 7 = c^ A{p, uj^)) / [mLOo) . At the peak, the signal-to- 
noise ratio, Eq. (|9|) divided by Eq. (^, is approximately 

2/2/^2. 
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